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t^Si' The effect of deformation on the two-neutrino double decay (2z//3/3-decay) for ground 

^ ■ state transition ''^Ge — > '^^Se is studied in the framework of the deformed QRPA with 

separable Gamow- Teller residual interaction. A new suppression mechanism of the 
CN I 2i^/3/3-decay matrix element based on the difference in deformations of the initial and 

final nuclei is included. An advantage of this suppression mechanism in comparison 
with that associated with ground state correlations is that it allows a simultaneous 
I description of the single f3 and the 2z//3/3-decay. By performing a detail calculation of 

^ ■ the 2z^/3/3-decay of '^^Ge, it is found that the states of intermediate nucleus lying in the 

. region of the Gamow- Teller resonance contribute significantly to the matrix element of 

00 I this process. 
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1 Introduction 



o 

^ ' There is a continuous attempt to improve the accuracy and reliability of the calculated 
^ ■ single beta {(3) and double beta decay nuclear matrix elements . Their values affect 
■ our understanding of astrophysical processes and the fundamental properties of neutrino, in 
^ . particular neutrino mixing and masses. 

The complexity of the calculation of the double beta decay matrix elements consists in 
the fact that it is a second order process, i.e., in addition to the initial and final nuclear states 
the knowledge of the complete set of the states of the intermediate nucleus is required. In 
solving of this problem different models as well as different nuclear structure scenarios were 
applied [1, 2, 3]. 

Double beta decay can occur in different modes: The neutrinoless mode (Oi//3/3-decay) re- 
quires violation of the lepton number and probes new physics scenarios beyond the Standard 
Model (SM) of particle physics [1, 2, 4, 5, 6, 7] . The two-neutrino decay (2z//5/3-decay) mode 
is process allowed in the SM, i.e., the half-life of this process is free of unknown parameters 
of the particle physics side [1, 2, 3, 4]. 

As the half-lifes of the 2i//5/5-decay have been already measured for about ten nuclei, 
the values of the associated nuclear matrix element can be extracted directly. It provides a 
cross-check on the reliability of the matrix element calculations. 

Since the nuclei undergoing double beta decay are open shell nuclei the proton-neutron 
Quasiparticle Random Phase Approximation (pn-QRPA) method has been the most em- 
ployed in the evaluation of the double beta decay matrix elements [1, 2]. In addition, this 
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approach has succeeded in reproducing the experimentally observed suppression of the 2u(3[3- 
decay transitions [10, 11, 12]. However, a strong sensitivity of the computed matrix elements 
on an increase of the strength of the particle-particle residual interaction in the 1"*" channel 
leads to a problem of fixing this parameter [1]. Thus various refinements of the original 
pn-QRPA have been advanced. 

In the last decade the extensions of the pn-QRPA, which include proton-neutron pairing 
(fuU-QRPA) [13], a partial restoration of the Pauli exclusion principle discarded in the pn- 
QRPA (renormalized QRPA) [14, 15] and higher order RPA corrections [16, 17], have received 
much attention in the context of the nuclear structure calculations for the double beta decay. 
In spite of the fact that these alterations of the pn-QRPA improved reliability of the evaluated 
matrix elements, further progress of the nuclear approaches dealing with double beta decay 
is needed. 

The deformation degrees of freedom of nuclei undergoing the 2z//3/3-decay were first con- 
sidered within Nilsson model with pairing [18]. Bogdan, Faessler, Petrovici and Holan cal- 
culated the 0z//9/9-decay in the pairing model including deformation and rotation [19]. The 
first QRPA calculation of the 2z//9/9-decay matrix elements in a deformed Nilsson-BCS basis 
were presented in Ref. [20]. The authors did not take into account particle-particle interac- 
tion of the nuclear Hamiltonian and assumed initial and final nuclei to be equally deformed. 
The effects of nuclear deformation were considered also within the SU{3) scheme [21], which 
has been found successfuU in describing the heavy rotational nuclei. The SU (3) scheme is 
a tractable shell model theory for deformed nuclei, which requires a severe truncation of 
the single particle basis. This approach was used for calculation of the double beta decay 
half-lifes of different heavy nuclear systems [22] and the predictions were found to be in good 
agreement with available experimental data for ^^^Nd and '^^^U [8]. The effect of the defor- 
mation of the nuclear shape on the two-neutrino double beta decay matrix element has been 
discussed in details within a method developed by Raduta, Faessler and Delion [23]. The 
authors used angular momentum projected single particle basis having the energies close to 
those of Nilsson levels. The Gamow- Teller states were generated with help of the spherical 
proton-neutron QRPA with a good angular momentum quantum number within the consid- 
ered basis. The results were presented for the 2i//3/3-decay of ^'^Se. It was shown that the 
deformation affects significantly the 2z//5/5-decay matrix element. The deformation effect on 
the double Gamow- Teller matrix element of ^^^Mo were investigated in the Hartree-Fock- 
Bogoliubov (HFB) framework in Ref. [24]. It was noticed that there is a necessity of an 
appropriate amount of deformation in the HFB intrinsic state to reproduce the experimental 
2i//3/3-decay half-life. 

There is an interest to study the effect of deformation on the double beta decay matrix 
elements within the deformed QRPA [25, 26], which allows unified description of the 2vf3j3- 
decay in spherical and deformed nuclei. This nuclear structure aproach has been found 
succeessful in description of the single /3-decay transitions of medium and heavy nuclei a 
long time ago. In the first applications only the particlc-liolc terms of the Gamow- Teller 
force were taken into account [25, 26]. It was supposed that the particle-particle terms have 
minor effect on the Gamow- Teller strength function. However, from the spherical QRPA 
calculations we know that the particle-particle force plays an important role for describing 
the and /3/3-processes [10, 11, 12]. Recently, the importance of the particle-particle 
interaction has been confirmed also in the deformed QRPA treatment of the Gamow-Teller 
strength distributions [27]. A strong sensitivity of the single /3-decay characteristics to the 



2 



nuclear shape, RPA ground state correlations and pairing correlations were manifested. 

The aim of the present paper is to investigate the effect of nuclear deformation on the 
double Gamow-Teller matrix element of the 2i//3/3-decay within the deformed QRPA. We 
present here as an example numerical results for the 2i//9/3-decay of ^^Ge, which is a prominent 
transition due to ongoing and planned /?/? experiments for this isotope [28, 29]. A subject of 
our interest is also the and /5+ distributions in ^^Ge and "^^Se, respectively. The problem 
of fixing of the nuclear structure parameters will be addressed. 

The paper is organized as follows. The basic elements of the deformed QRPA formahsm 
are reviewed in Section 2. In Section 3 we establish our choice for deformation, pairing and 
force parameters. Then the results for the 2z//?/?-decay ground state transition '^^Ge — Se 
are analyzed in context of the Gamow-Teller strength distributions of ^^Ge and '^^Se. The 
final conclusions and remarks are pointed out in Section IV. 

2 Formalism of the deformed QRPA 

The theoretical approach is based on the deformed proton-neutron Quasiparticle Random 
Phase Approximation with separable proton-neutron residual interaction, which is relevant 
for the allowed Gamow-Teller transitions [25, 26, 27]. This approach allows for an unified 
description of the 2z//3/3-decay process in deformed and spherical nuclei. 
The total nuclear Hamiltonian takes the form 

H ^ Ho + Hint- (1) 

Hq denotes the Hamiltonian for the quasiparticle mean field described by a deformed axially- 
symmetric Woods-Saxon potential [36] 

Hq^YI Erolpr^rpr {-r ^P,n), (2) 

TP 

where Er are the quasiparticle energies. a+^^ (o-rp^) is the quasiparticle creation (annihilation) 
operator. The p (n) index denotes proton (neutron) quasiparticle states with projection Qp 
(fin) of the full angular momentum on the nuclear symmetry axis and parity TTp (7r„). The 
index p (p = ihl) represents the sign of the angular momentum projection Q. We note that 
intrinsic states are twofold degenerate. The states with Q^- and — have the same energy 
as consequence of the time reversal invariance. We shall use notation such that is taken 
to be positive for states and negative for time reversed states. 

The method includes pairing between like nucleons in the BCS approximation with fixed 
gap parameters for protons, Ap, and neutrons, A„. 

= ^(£, - A,)2 + A2, (3) 

where Sr are single-particle energies. A,- is the gap, is the Fermi energy. The Bogolyubov 
transformation, which defines the quasiparticle representation, is given by 




Here, c^^^ and c^-p^ are creation and annihilation operators for particles, respectively. ~ 
indicates the single nucleon state obtained by the time reversal operator T — Ke''^"'^^^'^ — 
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(l 0^) '^^^^ state is given in the Appendix A [37]. The occupation probabihties v. 
are written as 



and the proton and neutron number equations are 

Z^2Y,vl N^2Y,vl (6) 

p n 

where Z and N are the numbers of protons and neutrons, respectively. 

The residual interaction part Hint of nuclear Hamiltonian in Eq. (1) contains two terms 
associated with particle-hole (ph) and particle-particle (pp) interaction: 

ft'=0,±l JS'=0,±1 

The operators /9~ and P~ are ph and pp components of the spin-isospin T"'"cr, namely, 

pPpUpn 

Pr ^ <PPp\^^^K\npn> cl^/npn, = {Pr)^ ■ (8) 

PPpTipn 

The ph and pp forces in Eq. (6) are defined to be repulsive and attractive (x, > 0), 
respectively, reflecting the general feature of the nucleon-nucleon interaction in the J'" = l"*" 
channel. The explicit form of the matrix element < ppp\T'^crK\npn > is presented in the 
Appendix A. 

After neglecting the scattering terms a^pp^Unp^ and CLnp^^pp-p quasiparticle representa- 
tion of Hint takes the form 

Hint = X E YXM^^K{j)+aK{i)aK{mA\i,K)AXj.K)+AQ,K)A{i,K)) 

K=0,±l ij 

+ {aK{t)(rK{j) +^K{tpK{jmA^{t,K)A{j,K) + A{j,K)A^{t,K))] 

H E[-(^^W7f;.(i)+^^(^)^^(i))(^^(^'^MHJ,^)+^(j,^M(^,^)) 

i<:=o,±i ij 

+ {nK{i)nK{j)+WK{i)WK{j)){A\i,K)A{j,K) + A{j,K)A\i,K^^^^ (9) 

with 

CrK{i) = < PPp\T^CrK\npn > UpVn, axii) =< PPp|T+(7x|np„ > VpUn, 

T^xii) = < PPp\T'^CrK\npn > UpUn, Wxii) =< PPp\T'^CrK\npn > VpVn- (10) 

AI and Ai are the two-quasiparticlc creation and annihilation operators 

A^i, K) = ala^^, A^{i, K) = aj^^at^^, K) = {A\i, K))l (11) 

The quasiparticle pairs i and i are defined by the selection rules flp — fln — K and fln — fip — 
K, respectively, and TTpTr^ = 1. 
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The above considered model Hamiltonian includes terms with K = 0, ±1 and describes 
J^K = 1+1, 1+0 excitations. In the laboratory frame the proton- neutron QRPA phonon 
wave functions for Gamow- Teller excitations in even-even nuclei have the form 



\lM{K),m> = 
\lM{K),m> = 




tl 



Vl,^{<l>,9,7P)Q^^\rpa> (X = 0), 



(12) 



where \rpa > denotes the QRPA ground state. The intrinsic states are generated by the 
phonon creation operator 



Qrni ^ ^[Xr^^A^t, K) - Y^^A{l, K)\. 



(13) 



In the case K = ±1 {K = 0) the sum in Eq. (13) includes all bound and quasibound two- 
quasiparticle spin-projection-flip (non-spin-projection-flip) configurations. We note that 
the K — —1 and K — 1 modes are related to each other through time reversal and are 
degenerate. 

The excitation energy Uk and the amplitudes and Y^]^ of the phonon Q]^^ are 

obtained by solving the RPA matrix equation 



fAiK) BiK)\ (X 
\B{K) A{K)) \Y, 



K 



1 

-1 



K 



K 



(14) 



where 



Aj(-f^) = ^Aj + 2x[(7;f (i)(7K(j) + (TK{i)'7K{j)\ - 2K[7r/f (i)7rK(j) + TT K{i)T^ K{j)\, 

Bij(K) = 2x[aK(i)aKij) + ax(i)(7x(i)] + 2KK(i)7f^(j) + 7f;f(i)7rK(j)] (15) 



with £i = Ep + En the two-quasiparticle excitation energy. 

An advantage of using the separable forces is that the RPA matrix equation reduces to 
a homogeneous system of only four equations for the four unknown norms N^r, N^, N^^ and 
N^, which is much easier to solve in comparison with the full diagonalization of RPA matrix 
of large dimension. The corresponding secular equation is given by 



1 + xi-p^l + Rf,: 



with 



K 



1717 J 



x(-Pj| + i?l, 
1 + x{-P^, + Rt) 

. pK _ TDK \ 



x{-Pl'. ~ RD 



<P^. + R, 

1 + k{P^^ - R 



kP. 



_|_ D/ 

K _ r)K_\ 
K _ T?K 

■KIT 



P-KTi) 



<P^. + Rt) 
<P.% + Rt) 
<P.% - RD 

1 + K{Pg - 



P^- = 2V 



K 



Si' 



a(i)a{i) 



(a, q; = (7, cr, tt, tt). 



K 



The forward and backward amplitudes are written as 



2N _ N- -'TT _ 

;7r^[ xi(^K{i) + (^Kii)^-) - «(7rx(i)-^ + 7fK(0^) 



K 



— 2Na- Ng- Nt, Ng- 



= 0, 
(16) 
(17) 



(18) 
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where for the norms N^^ N^, and N^^ we have 

j 3 

The normahzation factor is determined from the condition 

< rpa\m,Q^^]\rpa >= Y^i^rKX^K ' y^l^yl^) = 1- (20) 

i 

We note that the QRPA equations are calculated separately for different values of K and 
that the solutions for X = +1 and K — —1 coincide to each other due to considered axial 
symmetry. 

The (3~ and (3^ transition amplitudes from 0^ even-even initial nuclear state to a one- 
phonon state in odd-odd final nucleus are expressed by 

<lM(X),m|/3^|0+,.> = Y>i^XrK + a{z)Y^l 

i 

<IM{K)M^PmKs> = E[^(^)^1^ + ^(^)^^S- (21) 

i 

Here, |0^^ > denotes the correlated RPA ground state in the laboratory frame. The 
transition operators in Eq. (21) are related with intrinsic operators in Eq. (7) as follows: 

/3^ = E^mm(<^'^'WJ- (22) 

For the (3 strength function we have 

BgtH = E El<l(^),^ll/3-||0+,. >|'5(a;-a;K), 

K=0,±1 m 

B^A^) - E El<l(i^),^ll/3+l|0+,. >r<^(u;-a;x). (23) 

K=0,±1 m 

Prom the amphtudes one obtains, straightforwardly, the total strengths: 

Sgt = E El<i(^)'^ll/^"llo.';.. 

K=0,±1 m 

s^T = E EI<i(^)'^ll^'"llo.':..>l' (24) 

K=0,±1 m 

and for the Ikeda sum rule we get 

Sgt-S^t = E [Mi)?-{aKm 

K=0,±1 

= E E \<PPp\r-^^K\npn>?{vl-vl)^?,{N-Z). (25) 

K=0,±1 PPpUpn 
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The factor 3 comes from the sum over K, i.e., the contribution from each component K is 
equal to (N-Z). In deriving the above expression we used the closure condition for QRPA 
states, assumed the completeness relation for single particle states 

^\TPr><TPr\ = l (26) 

TpT 

and applied Eq. (6). We note that if a truncated single particle basis is considered in respect 
to that the full single particle basis of the Woods-Saxon potential, the condition in Eq. (26) 
is violated and as a consequence the Ikcda sum rule as well. 

The inverse half- life of the 2i//3/3-decay can be expressed as a product of an accurately 
known phase-space factor and the Gamow- Teller transition matrix element in 
second order: 

- 0.":..)]"' = i9AY \M'Jt\'- (27) 
The contribution from the two successive Fermi transitions is safely neglected as they come 
from isospin mixing effect [30]. Within the deformed QRPA approach the double Gamow- 
Teller matrix element M^.^ for ground state to ground state 2i//5/3-decay transition takes the 
form 

m2u _ X- \- < 0/ II r II HK),mf >< l{K),mf\liK),m, >< l{K),m, \\ r || 0+ > 

rriimf K=0,±1 ^ K ) I 

(28) 

The sum extends over all l"*" states of the intermediate nucleus. The index i (f) indicates that 
the quasiparticles and the excited states of the nucleus are defined with respect to the initial 
(final) nuclear ground state |0^ > (|0y >). The overlap is necessary since these intermediate 
states are not orthogonal to each other. The two sets of intermediate nuclear states generated 
from the initial and final ground states are not identical within the considered approximation 
scheme. Therefore the overlap factor of these states < l{K),mf\l{K),mi > is introduced in 
the theory [31]. Its takes the form 

< l{K),mf\l{K),m, >= iK^Kk - KkYi?^] T^hU < BCSf\BCS, > . (29) 

kif 

The factor TZiji^, which includes the overlaps of single particle wave functions of the ini- 
tial and final nuclei is given explicitely in the Appendix B. There, a detailed derivation of 
< l{K),mf\l{K),mi > together with the overlap factor of the initial and final BGS vacua 
can be find as well. In the spherical limit the value of the BGS overlap factor is about 0.8 and 
it is worth mentioning that it was comonnly neglected in the double beta decay calculations 
[11, 12, 33, 34]. 

We note that we neglected the overlap matrix elements between the intermediate states 
generated from initial and final nuclei with different K as the corresponding factors are 
very small due to different structure of corresponding RPA configurations. Thus the spin- 
projection-flip {K = ±1) and non-spin-pro jection-fiip [K — 0) excitations contribute co- 
herently to the 2i//3/3-decay matrix element Mq^. 

3 Calculation and Discussion 

The formahsm described in Section II is used for the calculation of the 2i//3/3-decay ground 
state transition "^^Ge — > '^^Se. The results are obtained with a deformed, axially symmetric 
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Woods-Saxon potential [35]. The deformation independent Woods-Saxon parameters (well 
depth, skin thickness, radius and spin-orbit constants) are taken from [36]. This parame- 
terization of the Woods-Saxon potential was used previously in different RPA calculations, 
where a good agreement was achieved with experimental data, in particular for single Ml 
transitions at low energy, observed in (e, e') and (7,7') experiments [38, 39]. 

The quadrupole-hexadecapole parameterization of the nuclear surface is considered: 

R(e) = Rn[l + ^2Y2o(0) + P4Y4o(e)l (30) 

where Rn is the normalizing constant {R{0 = 0) = 1) and Ylq (L = 2,4) are spherical 
harmonics. The spherical limit is achieved for P2 = Pa = 0. 

The deformation parameter /?2 can be deduced from nuclear electric quadrupole labo- 
ratory moment {Qp — —{7/2)Q, Q and Qp are laboratory and intrinsic quadrupole mo- 
ments, respectively) or extracted from values based on measured E2 probability {Qp = 

Y^T67rB(£^2)/5e^, the sign can not be extracted) via the intrinsic quadrupole moments Qp-. 

{tc is the charge root mean square radius). However, the available 
experimental data associated with the above two approaches [43, 44, 45] lead to a different 
quadrupole deformation both for ^^Ge and "^^Se. In particular, the quadrupole moments 
measured by Coulomb excitation reorientation method [44] , which yields also the sign of the 
intrinsic quadrupole moment, imply P2 — 0.1 (^^Ge) and P2 — 0.16 (^^-S'e) [41], but from 
the measured values of B{E2) strengths [45] one finds \(32\ = 0.26 (^^Ge) and |/92| = 0.29 
(^^Se) [41]. A lack of accurate experimental information on the deformation of "^^Ge and 
'^^Se suggests that it is necessary to study the associated 2z//3/3-decay matrix element as a 
function of deformation parameters of both initial and final nuclei. 

We note that there are theoretical calculations of deformation parameters performed 
within relativistic mean field theory [40], which imply 

76(2-6 : p2 = 0.157, /34 = 0.019, '^^Se : p2 = -0.244, (3^ = -0.028. (31) 

These results are in good agreement with predictions of the macroscopic-microscopic model 
of Moller, Nix, Myers and Swiatecki [42] as well as with nonrelativistic Hartree-Fock calcu- 
lations [41]. 

In the deformed QRPA calculation we assume a truncated model space by considering 
only single particle states with maximal allowed value of the asymptotic quantum number 

equal to 5. In contrast to many other microscopic calculations, energies of single particle 
levels are not shifted but taken exactly as provided by the deformed Woods-Saxon potential. 

The BCS equations are solved for protons and neutrons. The proton and neutron pairing 
gaps are determined phenomenologically to reproduce the odd-even mass differences through 
a symmetric five-term formula [46]. For both spherical and deformed shapes of the studied 
A=76 nuclei we consider the same values of the gap parameters, which are given by 

^^Ge : Ap = 1.561 MeV, A„ = 1.535 MeV, ^^Se : = 1.751 MeV, A„ = 1.710 MeV. 

(32) 

In Fig. 1 the calculated proton and neutron occupation probabilities close to the Fermi 
level are presented for '^^Se. We see that both for protons and neutrons the BCS solutions 
associated with spherical and deformed nuclear shapes differ significantly each from other. 
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This deformation dependence of the BCS results is expected to have strong impact on the 
QRPA solution. 

The calculation of the QRPA energies and wave functions requires the knowledge of the 
particle-hole x particle-particle k strengths of the residual interaction. The optimal 
value of X is determined by reproducing the systematics of the empirical position of the 
Gamow- Teller giant resonance in odd-odd intermediate nucleus as obtained from the {p, n) 
reactions [47, 48]. The parameter n can be determined by exploiting the systematics of single 
/3-decay feeding the initial and final nuclei. In Ref. [47] the strengths of the particle-hole and 
particle-particle terms of the separable Gamow-Teller force were fixed as smooth functions 
of mass number A by reproducing the /?-decay properties of nuclei up to A=150 within the 
spherical pn-QRPA model. The recommended values of x ^^'^ of R-ef. [47] are as follows: 

X = 5.2/^°-^ MeV, K = 0.58/^°-^ MeV. (33) 

In the presented 2z//3/3-decay study of the ground state transition ^^Ge — '^^Se we use the 
recommended value for x {x — 0-25 MeV) and k is considered as free variable. 

The QRPA calculations for the = 0, ±1 states are performed by following the procedure 
described in Section 11. The number of zeros of the non-linear RPA secular equation in Eq. 
(16) is equal to the number of the two-quasiparticle configurations included in the microscopic 
sums (17). In the case of = [K = ±1) the dimension of the configuration space is of 
the order of 450 (900) quasiparticle pairs excitations. To each zero of the RPA secular 
equation corresponds one RPA energy o;^. It is worthwhile to notice that by introducing 
the deformation degrees of freedom the zeros are not found in each subinterval of energy 
uj determined by the two subsequent two-quasiparticle excitation energies and that in some 
subintervals two or more zeros are found. This situation is illustrated in Fig. 2. The vertical 
axis is the left-hand side (l.f.s.) of Eq.(16) and the horizontal one is the energy ijJk=i in 
the range 6 MeV < uji.-i < 8 MeV . The upper (a) and lower (b) subfigures show results 
for A=76 obtained within the spherical and deformed QRPA, respectively. We note that a 
similar phenomenon was found also in Ref. [49] in different context. 

The main drawback of the QRPA is the overestimation of the ground state correlations 
leading to a collapse of the QRPA ground state, which may lead to ambiguous determination 
of the 13 and /3/3-decay matrix elements. The QRPA collapses, because it generates too many 
correlations with increasing strength of the attractive proton-neutron interaction as a result 
of the quasiboson approximation used. An interesting issue is the dependence of the position 
of the collapse of the QRPA solution on the deformation degrees of freedom. In Fig. 3 the 
energies of the lowest K — {) and = 1 states in "^^As calculated within the deformed QRPA 
in respect to the "^^Se ground state are plotted as a function of particle-particle interaction 
strength k. The calculations were performed for different values of quadrupole deformation 
parameter (32-, namely (32 = —0.25, —0.10, [32 = 0.0 and (32 = 0.10,0.25 representing oblate, 
spherical and prolate shapes of ^^5*6, respectively. We note that in the spherical limit the 
energies oi K — and = 1 states of the intermediate nucleus coincide to each other. Fig. 
3 shows that for a deformed nucleus the collapse of the QRPA solution appears for a shghtly 
smaller value of k. In the case of strongly deformed '^^Se this effect is more pronounced in 
comparison with the less deformed one. Thus the deformation degrees of freedom do not 
improve the stability of the QRPA solution. The violation of the Pauli exclusion principle 
affects the QRPA results strongly in the range of particle-particle interaction k above 0.07 
MeV. 
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The half-life of the 2z//3/3-decay of ^^Ge is known with high accuracy from the Heidelberg- 
Moscow experiment, in particular T^J^ = [1.55 ± 0.01{stat)to,llsyst.)] x 10^^ years [28]. All 
existing positive results on the 2i//3/3-decay were analyzed by Barabash [9], who suggested to 
consider for the ground state transition ^^Ge — > '^^Se the average value T^J^ — l-43lo!o7 ^ 10^^ 
years. By using of Eq. (27) and the knowledge of the kinematical factor G^'^ [G'^^(^^Ge) — 
1.49 X 10~^° year~^ MeV"^] from the 2z//5/3-decay half-life of ''^Ge one can deduce the absolute 
value of the nuclear matrix element {M^rf'^^^'l equal to 0.138 MeV^-^ by assuming = 1.25. 
If the value of the axial coupling constant qa is considered to be unity, the value of {M^rf'^^^' \ 
deduced from the average half-hfe for A=76 is larger, 0.216 MeV~^. 

In Fig. 4 the effect of the deformation on the 2i//3/3-decay matrix element Mq^j. is analyzed. 
The results for the 2z//9/9-decay matrix element are displayed as function of the particle- 
particle strength k. The curves drawn in subfigures a), b) and c) [subfigures d), e) and f)] 
correspond to the case both initial and final nuclei are oblate [prolate]. The two horizontal 
hues represent """^^ = 0.138 MeV'^ (qa = 1-25) and Maf^''^' = 0.216 MeV'^ (qa = 
1.0). By glancing Fig. 4 we note that within the whole range of k there is only a minimal 
difference between values of corresponding to the same value of A/92, which is defined 
as 

AP2 = \P2C'Ge)-P2C'Se)\. (34) 

In addition, one finds that by increasing the value of A/32 the suppression of M^.^ becomes 
stronger within the range MeV < n < 0.06 MeV. This a new suppression mechanism of 
the 2i//3/3-decay matrix element namely, depends strongly on difference in deformations 
of parent and daugter nuclei. 

One might ask what is the origin of this suppression. In Fig. 5 this point is clarified by 
presenting the overlap factor of two BCS vacua for different values of /32(™Ge) as function 
of the quadrupole deformation of '^^Se. We see that for a given /?2(^^Ge) the curve has a 
maximum for /32(^^Se) — /32(^^Ge) and with increasing difference in deformations of initial 
and final nuclei, i.e., A/92, the value of the BCS overlap factor decreases rapidly. From the Fig. 
5 it follows that oblate-prolate (or prolate-oblate) 2u f3 f3-decaj transitions are disfavoured 
in comparison with prolate-prolate or oblate-oblate ones for the same absolute values of 
associated quadrupole deformation parameters. In particular, by assuming the P2 and P4 
parameters of Ref. [40] we obtain < BCS(^^Ge)\BCS(^^Se) >= 0.009. This suppression is 
too strong to achieve agreement with experimental data for the calculated M^j,. However, by 
changing the deformation of '^^Se from oblate to prolate {(32C^Se) = 0.244) we end up with 
significantly larger value of the BCS overlap factor: < BCSC^Ge)\BCSC^Se) >= 0.61. 

The presented suppression mechanism is expecting to work also for other 2z//3/5-decay 
transitions. In Fig. 6 we present BCS overlap factor for the 2i//9/9-decay of ^^Ge, ^^°Mo, 
^^^Te and ^^^Xe. The quadrupole deformation of the initial nucleus was taken to be (32 — 0.1 
and (32 related with the final nucleus was considered to be a free parameter within the range 
—0.4 < (32 < 0.4. The pairing gaps, which entered the BCS equations, are given in Eq. (32) 
for A=76 and for A=100, 130 and 136 are given by 

^°°Mo : Ap = 1.612 MeV, A„ = 1.358 MeV, : Ap = 1.548 MeV, A„ = 1.296 MeV, 

^^°Te : Ap = 1.043 MeV, A„ = 1.180 MeV, ^^°Xe : Ap = 1.299 MeV, A„ = 1.243 MeV, 
^^^Xe : Ap = 0.971 MeV, A„ = 1.408 MeV, ^^^Ba : Ap = 1.245 MeV, A„ = 1.032 MeV{35) 

We see that the behaviour of the overlap factor for all considered nuclear systems is quahta- 
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tively the same. The maximum of this factor appears for equal quadrupole deformation of 
the initial and final nuclei. With increasing value of \l32{initial) — I32{final)\ the value of the 
overlap factor strongly decreases. This fact implies that the deformation of nuclei plays an 
important role in the calculation of the 2i//9/3-decay transitions and should be known with 
high reliability. 

An interesting issue is whether it is possible to obtain information about the deformations 
of ^^Ge and '^^Se by studying the /3~ and [3^ strength distributions in these nuclei, respec- 
tively. For different values of quadrupole parameter (^^ they are presented in Figs. 7 and 8. 
In order to facilitate comparison among various calculations the Gamow-Teller distributions 
are smoothed with a Gaussian of width 1 MeV. The assumed strengths of the particle- hole 
and particle-particle interaction strength were x = 0.25 MeV and k, = 0.028 MeV [47]. By 
glancing at Fig. 7 we see that the position of the Gamow-Teller resonance is not sensitive to 
effect of deformation and that only a httle strength is concentrated above the region of the 
Gamow-Teller resonance. However, there are significant differences between the spherical 
and deformed strength distributions in general [41] especially in the case of prolate shape of 
'^^Ge (see Fig. 7a), which is favoured by experimental data [44]. The thick solid line in Fig. 
7 corresponds to data obtained from charge exchange ^^Ge(p, nY^As reaction at 134.4 MeV 
[50], which were folded with 1 MeV width Gaussians as it has been done for the theoretical 
results. We see that the presented calculations are in a very good agreement with data for 
the position of Gamow-Teller resonance. It confirms the correct choice of the particle-hole 
strength of the nuclear Hamiltonian {x = 0.25 MeV) in our calculation [47]. From Fig. 7 
it also follows that the strength in the peak of the experimental Bqj, distribution is better 
reproduced by the QRPA calculations with only sligtly deformed mean field. 

The calculated strength distributions in "^^Se are displayed in Fig. 8. They have been 
folded with F = 1 MeV width Gaussians. We see that the effect of deformation on the f3~^ 
distribution strength for Gamow-Teller transition operator is more apparent as it is in the 
case strength distribution in ^^Ge. In Fig. 8 wc present also the experimental data, which 
are represented by solid points. They were obtained from charge exchange ^^5'e(n,p)'''M6" 
reaction [51]. We note that the experimentally extracted strength distribution with the help 
of two plausible methods differ from each other by more than a factor of 2 [51]. Hence the 
experimental data must be considered to be of a qualitative nature only. 

There arc two suppression mechanisms of the 2z//9/9-decay matrix element. Already long 
ago Vogel and Zirnbauer showed that Mq^ is strongly suppressed when a reasonable amount 
of particle-particle interaction is taken into account [10] , what is achieved close to a collapse 
of the QRPA solution. In this paper we show now that the 2z//9/3-decay matrix element Mq^ 
depends strongly on the difference in deformations of the initial and final nuclei. A possible 
criterium to decide, which deformation for the initial and the final nucleus in the double 
beta decay should be used could be the requirement for a common description of both single 
/3 and 2i//3/3-decay within the same nuclear Hamiltonian. In order to clarify this point three 
different cases are considered: 

i) Spherical shapes of both initial and final nuclei, i.e., f32{^^Ge) = and f32C^Se) = 0. 

ii) Deformations of parent and daughter nuclei: f32C^Ge) = 0.1 and f32C^Se) = 0.266. 

iii) Quadrupole deformations of the ground state of ^^Ge and '^^Se: (32C^Ge) =0.1 and 
P2C^Se) = 0.216. 

The choice of the deformation parameters in cases i) and ii) do not contradict the data 
[43, 44, 45]. 
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In Fig. 9 the 2i^/?/?-dccay of ^^Ge is presented as function of particle-particle interaction 
strength k, for the above three cases. In the calculation determined by case i) (spherical 
nuclei) M^^^""'^- equal to 0.138 MeV'^ {qa = 1-25) is reproduced for k = 0.06 MeV. The 
corresponding set of nuclear structure input parameters in this calculation is denoted by 
capital letter A (sec Table 1). This value of k is significantly larger as the average value 
Kf3 = 0.028 MeV deduced from the systematic study of the single /9^-dccays [47]. One can 
hardly suppose that it is because ^^Ge and '^^Se are specific nuclear systems. They are no 
indications in favor of it. It seems to be a general problem within the spherical QRPA: 
the calculated nuclear matrix elements reproduce the experimental 2z//3/3-decay half-lifes as 
a rule for particle-particle strength close to its critical value given by the QRPA collapse. 
In the case the calculation determined by case ii) the situation is different. M^rf'^^^' is 
reproduced for considerable smaller value of particle-particle strength k = 0.028 MeV, 
which corresponds to that of Homa et al. [47]. The letter B donotes the coresponding set 
of nuclear parameters, which is listed in Table 1. If we consider the effective axial coupling 
constant qa = 1, used often in the nuclei, M^"^^^^' is larger, namely 0.216 MeV~^. Then, a 
smaller difference between deformations of nuclei entering 2z//3/3-decay process is needed to 
reproduce this value for n ~ Kp. This is achieved with calculation of Mqt performed within 
the case iii) for k = 0.028 MeV (nuclear parameter set C). Nevertheless, one should keep 
in mind that the coupling strength of [47] has been adjusted by using of a deformed Nilsson 
single-particle model and qa = 1.25. Thus it is not necessarily the best possible choice. It is 
supposed that this prediction for k is more accurate for nuclei with shorter /5-decay half-lifes. 
We note that that the curves corresponding to cases ii) and iii) do change only weakly in 
the range < k < Kp. It means that the role of the particle-particle interaction within this 
interval is negligible. 

It is worthwhile to notice that for a large value of n of 0.07 MeV to 0.075 the agreement 
with [M^^^^^l deduced from the 2z//3/3-decay half-life is achieved as well and that in this 
case the sign of Mq^j, is negative. However, for negative value of the corrcspondance 
with from systematic studies of the single beta decay [47] is not achieved. Thus negative 
values of M^,^ are disfavored due to a complete disagreement with single /3-decay study of 
Homma et al. [47]. 

In what follows we shall discuss the single j3 decay characteristics corresponding to three 
different results in Fig. 4, denoted by letters A, B and C, for which the experimental half- 
life of the 2z//3/3-decay of ^^Ge is reproduced. The corresponding mean field and residual 
interaction parameters are given in Table 1. 

The (3~ and amplitudes calculated from the ground state of '^^Ge and "^^Se, respec- 
tively, associated with the K = {K = ±1) intermediate states of ™As are drawn in Fig. 10 
(Fig. 11). By glancing these figures we see that deformation (B and C sets of parameters) 
contribute to a fragmentation of the /3-decay amplitudes and that values of amplitudes 
associated with higher lying states of the intermediate nucleus (up to 20 MeV) are signifi- 
cant. It is worthwhile to note that the absolute values of = and K = ±1 amplitudes are 
equal each to other in case of the QRPA calculation with spherical mean field. The phases 
of these amplitudes are not of physical meaning. However, the relative phases of correspond- 
ing (3~^ and amplitudes plays a crucial role due to partial destructive character of the 
contributions from different intermediate states to M^,^. In the calculation these phases are 
fixed in such way that diagonal elements of the overlap matrix in Eq. (29) are taken to be 
positive. We note that the K — and K — ±1 contributions to M^,^ differ only slightly. 



12 



It is worthwhile to notice that a lot of (3^ strength in ^^5*6 is concentrated in the high 
energy region above 13 MeV (see Fig. 8), which can not be described within the shell model. 
There is a longstanding question whether in the calculation of the 2i//3/3-decay matrix element 
M'o'rp the transitions to higher lying states of the intermediate nucleus play an important role. 
There is so called "Single State Dominance Hypothesis" (SSDH), which assumes that only 
the lowest 1+ state of the intermediate nucleus is of major importance in the evaluation 
of M'o'rp. The SSDH can be realized, e.g., through cancellation among the higher lying 1"*" 
states of the intermediate nucleus. In Ref. [52] a possibility to decide experimentally about 
validity of the SSDH by measuring the single electron spectra and the angular distribution 
of the emitted electrons was discussed. This treatment of the SSDH is suitable for the 
2v j3 j3-dcc&y transitions with low-lying 1+ ground state of the intermediate nucleus. But this 
is not the case of the 2i//3/3-decay of "^^Ge. In Fig. 12 the QRPA model dependent study 
of this problem is presented by drawing the running sum of matrix element M'qj~{E) as a 
function of the 1+ excitation energy E^x of the intermediate nucleus. We see that in the 
case of the spherical calculation the main contribution to M'qj, comes from the states of the 
intermediate nucleus lying below 5 MeV scale and that there is a partial cancellation among 
contributions to M^.^ from higher lying states. In the case of deformed QRPA calculations 
there is a different situation. We see that important contributions to Mg,^(£') arise even 
from relatively high energies of the intermediate states around the position of Gamow-Teller 
resonance (about 12 MeV), which can not be described for a given nuclear system in the 
framework of the shell model. These results contradict the SSDH and indicate that a strong 
truncation in consideration of the complete set of the states of the intermediate nucleus is 
not appropriate. 



4 Summary and Conclusions 

For the 2i//3/3-decay ground state transition "^^Ge — > "^^Se a detailed study of nuclear matrix 
elements within the deformed QRPA with separable spin-isospin interactions in the particle- 
hole and particle-particle channels was performed. A new mechanism of the suppression of 
the double beta decay nuclear matrix elements was found. We showed that the effect of 
deformation on the 2i//3/3-decay matrix element is large for a significant difference in defor- 
mations of the parent and daughter nuclei and is not related with the increasing amount of 
the ground state correlations close to a collapse of the QRPA solution. The origin of this 
new mechanism is a strong sensitivity of the overlap of the initial and final BCS vacua to 
the deformations of the initial and final nuclei. This enters directly into the the overlap of 
the intermediate nuclear states generated from the initial and final nuclei via QRPA diago- 
nalization. A study of other double beta decay transitions as a function of this deformation 
dependent overlap indicates the general importance of the new suppression mechanism for 
the 2z//9/9-decay and 0z//?/9-decay matrix elements. 

Next, we showed that by assuming spherical nuclei there is a strong disagreement be- 
tween the particle-particle strength k needed to reproduce the 2z//3/3-decay half-hfe of ^^Ge 
and the k deduced in Ref. [47] from a systematic study of the single /3-decay transitions. By 
the new suppression mechanism of the 2z//3/3-decay matrix clement we showed that for defor- 
mations of the parent and daughter nuclei, which do not contradict the experimental data, 
a simultaneous description of both single /5 and 2i//3/3-decay is possible. We also discussed 



13 



the role of the axial coupling constant g^- 

We found that a significant part of the total /5+ strength is concentrated in the energy 
region above 10 MeV. A detail study of contributions to from different 1+ intermediate 
nuclear states showed that for the A=76 nuclei it is necessary to consider all l"*" states of 
up to energy 15 — 20 MeV. This fact disfavors theoretical studies of the 2i//3/3-decay of 
medium and heavy nuclei in models with a strongly truncated basis, like the shell model. 

The Gamow- Teller strengths and double beta decay matrix elements in medium and 
heavy nuclei continue to be a challenge for nuclear structure models. Further theoretical 
studies are needed to clarify the role of deformation in calculation of other double beta 
decay transitions, especially those including heavier nuclear systems, which are known to be 
strongly deformed. Future works should investigate the 2z//3/3-decay matrix elements also 
within different possible extensions of the deformed QRPA, e.g., those including proton- 
neutron pairing [55] . In this way one expects to develop a reliable many-body approach with 
well defined nuclear structure parameters for calculation of the Oi^/3/3-decay matrix elements. 
Their accurate values are highly required to determine the neutrino mixing pattern and to 
answer the question [53] about the dominant mechanism of the 0i//9/9-decay [54]. 

5 Appendix A 

The single-particle states are calculated by solving the Schrodinger equation with the de- 
formed axially symmetric Woods-Saxon potential, which parameterization is given in Ref. 
[36]. They are characterized by their energy Et, parity tt,- and by the projection Jl^ (r = p, n) 
of the full angular momentum on the nuclear symmetry axis. We use notation \ppp > and 
\npn > for protons and neutrons, respectively. |rp^ > represents proton (r = p) or neutron 
(r = n) state with quantum numbers fir and tIt- Pt is the sign of the angular momentum 
projection Q (p^ = ±1). We note that intrinsic states are twofold degenerate. The states 
with fi^ and — fir have the same energy as consequence of the time reversal invar iance. We 
shall use notation such that Pr is taken to be positive for states and negative for time reversal 
states. 

In order to solve the Schrodinger equation the eigenfunctions of a deformed symmetric 
harmonic oscillator are used as a basis for the diagonalization of the mean-field Hamilto- 
nian [35]. These states are completely determined by a principal set of quantum numbers 
{N,nz, A,fl), where N = Ux + Uz, = 2n.r + |A| and $7 = A + S. — 1 and — 1 
are number of nodes of the basis functions in the ^-direction and r-direction, respectively. 
A and S are the projections of the orbital and spin angular momentum on the symmetry 
axis z. The explicit form of single-particle harmonic oscillator wave functions in cylindrical 
coordinates (r, z, 0) can be found, e.g., in Ref. [37]. For a given shape of the nuclear surface, 
the shape of the deformed harmonic oscillator is automatically chosen in a way suitable for 
obtaining good accuracy with a smallest number of basis functions [35]. The deformation 
dependent cutt-off is chosen in such way as to assure numerical stability of the results. In 
our calculation we use 11 major shells. 

In cylindrical coordinates the eigenfunctions of states and time-reversed states in de- 
formed Woods-Saxon potential are expressed as follows: 

|tp^ = +1> = ^[6t2^^jAr,n„A„i7, = A^ + 1/2 > + 
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bhiajN, n,, A, + 1, = A, + 1 - 1/2 >] (36) 

and 

\fpr = +l> = \TPr = -l>^Y.[b'}^LnrW, n„-Ar,fl= -Ar- 1/2 > - 

b^NlnjN, n„ -A, - 1, = -A, - 1 + 1/2 >] (37) 

with A > 0. ~ indicates time reversal states. 

The single-particle matrix elements of the T'^ax operator are given by 

<pPp\T+aK=o\npn> = fepn„PpI][&S\l2.n^^Svia„ -&Sviap&!vinJ, (38) 

<pPp\r+aK=+i\npn> = -V^<^npf2„+i $1 ^'S"inp^'Svin„ for Pp^Pn^+1, 

= +V2(5n^n„+i J2 ^Svin^^'Svion for Pp = p„ = -1, 

= -^WM-\T.bNln/Nln^ /or = +1, p„ = -X,39) 

<pPp\r+aK=-i\npn> = v^f^n^nn-i ^!vinp^Svin„ M Pp = Pn = +1, 

Nuz 

= -V^5Qpf^„-l ^ &S\l'inp^SviQ„ /O?^ Pp^Pn^ -1, 

= E 41f2,&^in„ /or Pp = +1, p„ = -1. (40) 

The overlap of the proton and neutron single particle states of the initial (^4, Z) and 
final {A, Z + 2) nuclei are calculated by assuming that the corresponding sets of basis wave 
functions do not differ significantly each from other. Then we have 

< PfPpf\PiPpi > = ^np^np^^[bNlznp,bNlznp. ~ bNlzQpfbNlznp.]^ 

Nnz 

<nfPnf\nipni> = fe„^n„, E[^Al2.n„,^ivin„. - ^Svin„ (41) 

Nuz 

The index i (f) denotes that proton and neutron single particle states are defined with respect 
to the initial (final) nucleus. 

6 Appendix B 

As a consequence of considered many-body approximations the two sets of intermediate 
nuclear states generated from the initial and final ground states are not identical within the 
QRPA theory. Thus it is necessary to introduce the overlap factor of these states, which can 
be expressed with the help of intrinsic phonon operators as follows: 

< l(X),m/|l(X),m, >=< rpa/l Q7 Vpai > . (42) 

Here, the index i (f) indicates that the excited states of the nucleus are defined with respect 
to the initial (final) nuclear ground state \rpai > {\rpaf >). 
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In order to evaluate < l{K),mf\l{K),mi > we express the phonon creation operator 
Q^'^ in terms of creation and anninhilation phonon operators associated with the final 
nucleus. We have 

Qk ~ {(^mimjQK ^ruimfQ k\ (43) 

The coefficients of expansion armmf a-nd brmmf will be determined below. By inserting Eq. 
(43) into Eq. (42) we get 

<l{K),mf\l{K),mi> = '^P^'fl^K^QK^^l'^P^^i X^rmm'f 

m'f 

+ < rpaf\QK^Q'^]^\rpai > brmm'f] 

am,m'f < BCSf\BCSi > . (44) 

Here, we neglected overlap matrix element between the final ground state and the two- 
phonon state generated from the initial nucleus, which is considered to be small and should 

be not related with the studied quantity. In addition, we approximated the overlap of initial 
and final RPA ground state with the BCS ones. Thus the overlap factor of the intermediate 
nuclear states generated from the initial and final ground states is proportional to the overlap 
of initial and final BCS vacua. 

We proceed with the calculation of armm'f and brmm' j coefficients. The quasiparticle 
creation and anninhilation operators (a*^*-'''^, a*^*-*) [(a*^'^^''', a*--^^)] of the initial [final] nucleus are 
connected with the particle creation and annihilation operators )] by the 

BCS transformation [see Eq. (4)]. In addition there is a unitary transformation between 
particle operators associated with initial and final nuclei 

E <rPr\T'pr'>c[;)l 

^ <Tp.|TV>c^i^, (45) 

The overlap factors of the single particle wave functions of the initial and final nuclei < 
tPt\t' p' ^, > is given explicitely in the Appendix A. The above mentioned transformations 
allow us, by using the quasiboson approximation, to rewrite the boson operators of the initial 
nucleus with the help of the boson operators of the final nucleus: 

I' 

A^^\lK) = Y.i^u'A(f\t,K)-Su'A^f^^{l',K)]. (46) 
I' 

Assuming the definition of the quasiparticle pairs operator in Eq. (11) the factors TZu, and 
Sw can be expressed as 

= < PPpIp'Pp' > (4'^4'^ + < nPn\n'pn' > {u^Sn' + 4*^4'^), 

Siv = < PPpIpV > {ufvf - ufv'^J,^) < npn\n'pn' > {u^M' - uM')- (47) 

It is worthwhile to notice that in the limit initial and final states are indentical TZw — 1 and 
Siv = 0. 



.(i)t 
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By inserting Eq. (46) into the expression for the phonon operator of the initial nucleus 
[see Eq. (13)] and by exploiting the relations 

A^'HlK) = UKKQ7 + y:]iQ7\ (48) 

we find [31] 

IV 

+ Y;;]iSinxr^-Y;r^SnYi^% (49) 

By neglecting the terms proportional to Sui due to their smallness we end up with the overlap 
factor of the intermediate nuclear states given in Eq. (29). 

The overlap factor of the initial and final BCS vacua can be written as product of proton 
and neutron BCS overlap factors for a given anular momentum projection quantum number 
Vt: 

< BCSf\BCSi > = < BCSf{p)\BCSi{p) >< BCSf{n)\BCSi{n) > 

= n < BCSf(np)\BCSi(np) > n < BCSf(nr^)\BCSi(nr^) > -(so) 

where 

< Bcsf{n)\Bcsm > = < I n (4'^ + n (4^ + i > ■ (51) 

k=l 1=1 

Nq is the number of single particle states with the same value of quantum number Q. The 
same model space for protons and neutrons is assumed. By a direct calculation of the above 
matrix element one finds 

< BCSf{n)\BCSi{n) >= JJ 4^^ 114^^ 

k=l 1=1 

Nq Nainii) Nn{ni) 

+ E e^-S {D^'\rn,-n,)) U u^^ H 4^ 

m\,n\=l k=l 1=1 

ATn Nn{mi,m.2) Nci{ni,n2) 

+ E v^viilvSv^{D^^im„m,;n„n^ H 4'^ U 

m\,m2,n\,n2=l k=l 1=1 

Nn Nn 

+ ••• + (i^(^")(l,2,...,iVn;l,2,..,7Vn)) 114'^ II^^- (52) 

k=l 1=1 

Here, Hfe^i'"^'™^'' means that index k runs the values from 1 to except the values k = 
mi and k = 1712 (1 < mi < and 1 < 1712 < iVn)- -D''(mi, m2, m^; ni, n2, n^) 
denotes the determinant of matrix of rank r constructed of elements of the unitary matrix 
of the transformation between the initial and final single particle states with row indices 
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mi,m2, ...,771^ and column indices ni,n2, ...,nr. It is worthwhile to notice that by replacing 
all determinants in Eq. (52) with unity, i.e. the matrix of the transformation between 
the single particles associated with both nuclei is just unity matrix, we obtain a compact 
expression [20] 

< BCSfmBCS,{n) >= n(4^^«? + (53) 

fe=l 

However, this approximation is not justified and can lead to a significant inaccuracy in the 
calculation of Mq^ especially if there is a strong difference in deformations of the intial and 
final nuclei. 
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Table 1: Three sets of nuclear structure input parameters (A, B, and C) for which the 
calculated nuclear matrix elements M^r^ reproduces the experimental 2i//3/3-decay half-hfe 
of ^^Ge. v.o.f. denotes overlap factor of the initial and final BCS vacua. 
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Figure 1: The proton (a) and neutron (b) occupation probabilities close to the Fermi level 
for ^^jS'e. The spherical (sph, triangle up) and deformed (dcf, closed circle) BCS results 
correspond to pairing gaps in Eq. (32). The vertical lines denote Fermi energy for protons 
{Xp) and neutrons (A„). 
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Figure 2: The l.h.s. of Eq. (16) D{uj) is plotted as function of the energy uj for ^^Ge. The 
calculation has been performed within spherical (a) and deformed (b) QRPA for x — 0.25 
and K = 0. In case (b) the K=0 results are presented. 
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Figure 3: The energy of the lowest QRPA state as function of particle-particle interaction 
strength k,. The K=0 and K=l deformed QRPA results for '^^Se are presented in subfigures 
(a) and (b), respectively. 
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K [MeV] K [MeV] 

Figure 4: 2z//3/?-decay matrix element of ^^Ge as function of particle-particle interaction 
strength k. In the subfigures (a), (b) and (c) [(d), (e) and (f)] the results corresponding to 
oblate [prolate] deformation of both initial and final nuclei are presented. Please note that if 
the deformation of the initial and final nuclei is comparable, there is only minimal difference 
between the calculated values of M'^j-. By increasing difference in deformations of parent and 
daughter nuclei the suppression of is increased in the range MeV < k < 0.6 MeV. 
The two dashed horizontal lines correspond to M^rf'^^^ = 0.216 MeV~^ [qa = 1-0) and 
Mcf''''' = 0.138 MeV-^ {qa = 1-25). 
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Figure 5: The overlap factor of the initial and final BCS vacua as function of the quadrupole 
deformation parameter of ''^S'e. The results are presented for spherical (/?2 = 0.0), oblate 
{/32 = -0.25, -0.10) and prolate {(32 = 0.10, 0.25) deformations of ^Ge. 
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Figure 6: The overlap factor of the initial and final BCS vacua as function of the quadrupole 
deformation parameter of the final nucleus for double beta decay of ^^Ge, ^°°Mo, ^'^"Te. 
and ^^^Xe. The deformation parameter of the initial nucleus is chosen to be I32{'i'nitial) =0.1. 
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Figure 7: Gamov-Teller strength distributions in ^^Ge as a function of the excitation 
energy of the daughter nucleus E^x- The distributions associated with prolate and oblate 
deformation of ^^Ge arc compared with that obtained for spherical shape {P2 = 0.0) in 
subfigures (a) and (b), respectively. The recommended values of x and k, of Ref. [47] were 
considered: x = 0.25 MeV and k = 0.028 MeV . Experimental data (thick solid line) are 
from Ref. [50]. 
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Figure 8: Gamov- Teller strength distributions Bqj^ in ^^5*6 as a function of the excitation 
energy of the daughter nucleus E^x- The distributions associated with prolate and oblate 
deformation of ^^Ge are compared with that obtained for spherical shape {P2 = 0.0) in 
subfigures (a) and (b), respectively. It is assumed x — 0-25 MeV and k, — 0.028 MeV [47]. 
The data represented by the solid points are from Ref. [51]. 
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Figure 9: 2i//5/?-decay matrix element of ^^Ge as function of particle-particle interaction 
strength k. The solid (sph.) line corresponds to spherical shape of initial and final nuclei. The 
dashed (def. 1) and dot-dashed (def. 11) are asociated with a set of deformation parameters 
{H^^Ge) = 0.1, P2C^Se) = 0.266) and {P2C^Ge) = 0.1, f32C^Se) = 0.216), respectively. 
The points indicated by letters A, B and C determine k for which the value of M^^"^^'^' 
deduced from the 2z//?/?-decay half-life of ''^Ge is obtained. For = 1-25 {qa = 1-0) one finds 
M^u-exp ^ Q_^3g j^^Y~i (]Vf2^-"^'P = 0.216 MeV-^) by assuming Tlj^{^^Ge) = 1.43 x lO^i 
years [9]. 
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Figure 10: The transition amplitudes between the ground state of ^^Ge (^^5*6) and 

K — intermediate states as function of the excitation energy Eex of '^^As. The (3~ 
results obtained for input parameter sets A, B, C (see Table 1) are presented in subfigures 
a) [d)], b) [e)] and c) [f)], respectively. Please observe, that the scale is different for the P~ 
and the Z^"*" strength distributions. 
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Figure 11: The same as in Fig. 10 for transitions to X = ±1 intermediate states. 
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Figure 12: Running sum of the 2i//3/3-decay matrix element as a function of the exci- 
tation energy Eex in '^^As. In the upper figure (a) results corresponding to input parameter 
sets A (solid line) and B (dashed line) are presented. In the lower figure (b) results ob- 
tained with input parameter set C (dot-dashed line) are drawn. The dotted horizontal line 
in subfigure (a) [subfigure (b)] donotes the value of matrix element M^^^^^ = 0.138 MeV~^ 
{M^f^^'P = 0.216 MeV^) deduced from the experimental 2z//3/3-decay half-hfe of ^^Ge by 
assuming qa — 1-25 {qa — 1-0). 
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